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A set of molecular models for 78 pure substances from prior work is taken as a basis for 
systematically studying vapor-liquid equilibria (VLE) in ternary systems. All 33 ternary 
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O ■ mixtures of these 78 components for which experimental VLE data is available are studied 



by molecular simulation. The mixture models are based on the modified Lorentz-Berthelot 



■ combining rule that contains one binary interaction parameter which was adjusted to a 
single experimental binary vapor pressure of each binary subsystem in prior work. No 
adjustment to ternary data is carried out. The predictions from the molecular models of 
the 33 ternary mixtures are compared to the available experimental data. In almost all 
cases, the molecular models give excellent predictions of the ternary mixture properties. 
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1 INTRODUCTION 



In previous work of our group, a set of molecular models has been developed for 78 real pure 
fluids using the dipolar or quadrupolar two-center Lennard- Jones (2CLJD and 2CLJQ) po- 
tential PQ[2|- This model type has been proposed more than three decades ago [3], however, 
it is far from being fully exploited. Polar 2CLJ models consider the basic molecular in- 
teractions repulsion and dispersive attraction and also feature anisotropy and polarity in 
a simple way. 78 small molecules consisting of up to nine atoms that belong to different 
classes of real fluids, including noble gases, alkanes, halogens and numerous refrigerants, 
were modeled using that approach [H [2]. For many of the 78 molecules, the polar 2CLJ 
model is only a crude assumption. E.g., the asymmetry of molecules is neglected and the 
polar interaction is always aligned along the molecular axis. Also the polarizability, which 
is often assumed to be a crucial molecular property for thermodynamics, is only implic- 
itly considered by Lennard- Jones interaction sites. Furthermore, the internal degrees of 
freedom are neglected as the polar 2CLJ models are rigid. 

After showing that these simple molecular models give excellent results for vapor-liquid 
equilibria (VLE) of both the pure components and their binary mixtures [H Ej H], the 
aim of the present work is to investigate on a broad basis whether such models are fully 
transferable to VLE of ternary systems. 

Based on the 78 pure substance models [HIS], the unlike energy parameter was adjusted 
in previous work [ll[5l[6l[7] to the experimental binary vapor pressure for 267 binary systems 
in order to very accurately describe their VLE. The direct transferability of such models to 
higher systems was also shown by our group [HI El [7] with VLE predictions of five ternary 
mixtures. This work is extended here to 33 ternary systems. 

A few publications on molecular simulation results for ternary VLE are available from 
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different authors: Carrero-Mantilla and Llano- Restrepo [S] (N2 + CH4 + C2H6), Potoff 
and Siepmann ^ (N2 + CO2 + propane), Kamath and Potoff [lOj (CH4 + H2S + CO2), 
Hansen et al. [H] (N2 + O2 + CO2), Liu and Beck [I2] (CH4 + CO2 + C2H6), Nath et al 
[13] (C2H4 + 1-hexene + polyethylene), Lfsal et al. [H] (isobutene + methanol + MTBE) 
and Van't Hof [15j (CH4 + CO2 + C2H6 and CH4 + CO2 + propane). However, each 
of these publications is restricted to one or two ternary mixtures only. Note that there 
are additional works on ternary VLE by simulation [161 [13 [HI [13 [20l [21] , but they deal 
exclusively with model systems. 

To avoid an arbitrary selection of the studied systems, a combinatorial brute force 
approach was taken here. Theoretically, out of the N = 78 components modeled in [H [2] 
A^(A^ — 1)(A^ — 2)/6 = 76 076 ternary mixtures can be formed, but of course, by far not all 
of these systems have been studied experimentally. To our knowledge, VLE were measured 
only for a subset of 33 out of the 76 076 ternary systems, corresponding to 0.043%. In 
the present work, all these 33 ternary mixtures were studied. This is the largest set of 
ternary systems that was used so far to probe the application of molecular modeling and 
simulation to ternary mixtures. 

It would have been attractive to investigate VLE of multi-component mixtures with 
more than three components too, however, to our knowledge no experimental data exists 
for any mixture consisting of four or more of those 78 components. 

The simulation results from the present work are compared to experimental data and 
in most cases to the Peng-Robinson equation of state (EOS) [22j which was applied in the 
same predictive way, i.e. it was adjusted to the same binary data that was also used to 
adjust the molecular mixture models, cf. [1]. Due to the fact that the Peng- Robinson EOS 
is widely known, it is not described here, for details see [1]. 
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2 EXPERIMENTAL DATABASE 



Experimental data were predominately retrieved using Dortmunder Datenbank (DDB) [23] . 
which collects all publicly available mixture VLE data sets, covering more than a century 
of experimental work. According to DDB, for a subset of 33 of the potential 76 076 ternary 
mixtures experimental VLE data is available. They stem from 30 publications [2^- [53]. 
These 33 ternary systems include 35 of the 78 pure components; cf. Table [1] for the full 
component list including their CAS RN number for proper identification. Please note that 
the ASHRAE nomenclature is preferred in the following due to its brevity, despite its 
deficiencies [51] . 

The studied 33 ternary systems are listed in Table [2] together with a reference to the 
experimental VLE data. Of those 33 ternary mixtures, five have been modeled in previous 
work of our group [3 [H [7] but the resulting VLE data were published only partly. 

It can be argued that these 33 systems, being just 0.043% of the full combinatorial 
sample, were selected by the experimentalists due their technical or scientific importance. 

3 PURE FLUID MODELS 

As explained above, 35 polar 2CLJ molecular models, taken from Pi's], were used here. 
A list of the pure fluids is given in Table [H These are two spherical non-polar (ICLJ) 
models for Ar and CH4, three spherical dipolar (Stockmayer or ICLJD) models for R30, 
R30B2 and R32, furthermore 17 elongated dipolar (2CLJD) models which include CO 
and numerous refrigerants, and finally 13 elongated quadrupolar (2CLJQ) models which 
include N2, O2, alkanes, refrigerants and CO2. 

A detailed description of the polar two-center Lennard- Jones pair potential is provided 
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in [1] and not repeated here. Polar 2CLJ models have four parameters: size a, energy e, 
elongation L and either dipolar momentum ^ or quadrupolar momentum Q. Stockmayer 
models have a vanishing elongation, while the non-polar spherical LJ models have only a 
and e. Model parameters were adjusted in [H [2] to experimental pure fluid VLE data using 
global correlations of critical temperature, saturated liquid density and vapor pressure as 
functions of these molecular parameters [5^ . These pure substance model parameters 
are also not repeated here. It should be noted that a wide range of polar momenta are 
covered by the 35 pure substance models regarded in this work. Starting from a non- 
existent polar momentum in case of Ar and CH4, it ranges to up to 3.7104 D for the 
dipolar R30B2 and up to 16.143 DA for the quadrupolar RlllO. 

The advantage of these molecular models is their simplicity, which reduces simulation 
time considerably, and their accuracy: typically, the relative deviations between simulation 
and experiment are below 1 % for the saturated liquid density, below 3 % for the vapor 
pressure, and below 3 % for the enthalpy of vaporization [H [2]. They also have shown 
to reliably predict Joule-Thomson inversion curves for pure fluids and mixtures [571 EH] , 
covering a wide range of state points, and also transport properties [59l EOl Ell ESI ES] . 

4 MOLECULAR MIXTURE MODELS 

On the basis of pairwise additive pure fluid potentials, molecular modeling of mixtures 
reduces to modeling the interactions between unlike molecules. Unlike interactions consist 
of two different types here. On the one hand there are the electrostatic interactions (dipole- 
dipole, dipole-quadrupole, and quadrupole-quadrupole). These interactions are treated in 
a physically straightforward way, simply using the laws of electrostatics [1] . 

Repulsion and dispersive attraction are other interaction types and are present between 
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all molecules. If a mixture A + B + C is modeled on the basis of Lennard- Jones potentials, 
the knowledge of three pairs of unlike Lennard- Jones parameters is required: (Tab, ^ab and 
c"AC; cac as well as ctbc, ^bc- For their determination, the broadly used Lorentz-Berthelot 
combining rule is a good starting point [61]. However, introducing a binary interaction 
parameter ^ to adjust the unlike energy parameter Cjj 

(Tij = ((Ti + (Tj)/2, (1) 

and 

eij = i^e'e], (2) 

allows almost always for an optimal representation of the binary fluid phase behavior 

For VLE, it was shown in ^1] that ^ can be adjusted to a single experimental binary 
vapor pressure. Values for ^ are given in |3] for 267 binary combinations. Note that 
the present 33 ternary systems comprise 65 different binary subsystems, whereof 62 were 
covered in [1]. The three exceptions are N2 + R14, R125 + R161 and R134a + R161. 
It was abstained here to adjust the binary interaction parameter for these three binary 
subsystems to ternary VLE data, thus ^ = 1 was specified instead. We refrained here from 
adjusting the binary interaction parameter kij of the Peng-Robinson EOS for these. 

5 RESULTS AND DISCUSSION 

To assess the predictive quality of the mixture models, ternary VLE were determined by 
molecular simulation predominantly at state points for which a direct comparison to ex- 
perimental data is possible. Simulation details are given in the Appendix. The Grand 
Equilibrium method [78] was used for the VLE simulations, where temperature and liquid 
composition are the independently specified thermodynamic variables, while vapor pres- 
sure, saturated vapor composition, saturated densities and enthalpy of vaporization are 
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determined. In most cases, simulation results are presented that match exactly with the 
experimental bubble line composition. However, if it was found that there is a significant 
mismatch for the resulting vapor pressure with respect to the experiment, the simula- 
tive bubble line composition was altered so that both data sets are almost in the same 
temperature-pressure plane. 

As experimental densities and enthalpies are rarely available in the literature, only 
vapor pressure and saturated vapor composition were used for this assessment. It should 
be noted that saturated vapor composition data is available for 26 of the investigated 33 
ternary systems, for the remaining seven systems, cf. Table [2], only bubble line data is 
available. 

The results are presented here in ternary plots at constant temperature and pressure, 
cf. Figures [1] to [231 covering 23 of the 33 studied ternary mixtures. For the remaining 
ten systems less than two full experimental VLE data points are available for one pair of 
temperature and pressure so that the results are not presented in figures here, while the 
numerical comparison to experimental data can be found in the supporting information. 
Full numerical VLE simulation data are given in the supporting information, which also 
contains saturated densities and heat of vaporization from simulation. 

For all predicted VLE properties, an estimate of the statistical uncertainty is provided 
in the supporting information. Due to the fact that the error bars are mostly within symbol 
size, they were omitted in the figures. 

The present assessment was made on the basis of the resulting composition on the 
saturation lines which can directly be seen in comparison with the experimental data in 
the phase diagrams of this work. Note that the simulated vapor pressure in general does 
not match exactly with the presented pressure, but it is usually very close to it. The 
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numerical data in the supporting information allows for a direct comparison of the vapor 
pressure. 

Not for all studied systems, the experimental data is sufficient to assess the topology of 
the saturation lines in the the isobaric-isothermal phase diagrams. Most of those, for which 
this is possible, show a simple topology where one straight bubble line and one straight 
corresponding dew line connect two binary subsystems, e.g. Ar + CH4 + N2, cf. Figure [T], 
N2 + CO2 + R12, cf. Figure [13 or CO2 + R142b + R152a, cf. Figure [HI However, three 
mixtures have markedly curved phase envelopes, i.e. CH4 + CO2 + C2H6, cf. Figures [9], 
N2 + CO2 + C2H6, cf. Figure [in and R13 + R14 + R23, cf. Figure [201 AH three have one 
azeotropic subsystem P], which however, cannot directly be seen from the figures for the 
ternary systems shown here. The phase envelope is also curved for the mixture CH4 + N2 
+ C2H6, cf. Figure [71 which has a ternary critical point. Finally, Figure [T7| for the mixture 
C2H2 + C2H4 + C2H6, shows two pairs of straight saturation lines which also result from 
the azeotropic behavior of the subsystem C2H2 + C2H6 [1]. 

The temperature range covered in the present study is quite large, i.e. from 112 K (Ar 
+ CH4 + N2, cf. Figure [ID to 358.5 K (RIO + RlllO + R1120, cf. Figure [IHD. The same 
holds for the pressure range, i.e. from 0.07 MPa (RIO + R20 + R30, cf. Figure [IS]) to 12.4 
MPa (N2 + O2 + CO2, cf. Figure [TTl) . For most mixtures, experimental data was available 
only for one pair of temperature and pressure, however, for 11 ternary systems either two 
(10) or three (1) pairs were simulated. Thereby, the largest temperature interval was 50 
K (N2 + CO2 + C2H6, cf. Figure [121 and supporting information) and the largest pressure 
interval was 7.23 MPa (N2 + O2 + CO2, cf. Figure [TTl and supporting information). 

In general, it can be stated that the agreement between simulation and experiment is 
very satisfactory. Both qualitatively and quantitatively, an excellent match was found in 
almost all cases. 
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Only for the mixture N2 + R13 + R14, the VLE could not be simulated at the conditions 
for which experimental data [38] is available. Particularly the temperature of 77.8 K, which 
is well below the pure substance triple temperature of both R13 (92 K) and R14 (89.5 K), 
posed a problem during simulation as it is in immediate vicinity to solidification. It should 
be pointed out that the employed molecular models [H [2] were neither optimized nor 
evaluated with respect to the triple line. 

Also results from Peng-Robinson EOS with adjusted binary parameter kij are shown. 
This model is known to be a good correlation tool, making it a workhorse in process 
engineering, it performs satisfactory in many cases too. Within the 23 examples presented 
here, three cases can be identified where significant deviations to the experimental data 
can be seen, i.e. N2 + O2 + CO2, cf. Figure M, RIO + R20 + R30, cf. Figure M and 
R140a + R141b + R142b, cf. Figure [231 Finally, for Ar + N2 + O2, cf. Figure [1 the 
Peng-Robinson EOS deviates substantially from the remaining two data sets. 

6 CONCLUSION 

It was shown that molecular modeling and simulation is a reliable and robust approach 
to obtaining VLE of ternary mixtures. To verify this, a total of 33 ternary mixtures were 
studied by molecular simulation with the Grand Equilibrium method. This method was 
found to be well suited for simulations of ternary VLE. 

The pure component models used in the present study were adjusted to pure component 
VLE data in previous work. For the binary mixtures, one state-independent parameter was 
adjusted to binary VLE data in other previous work as well. Due to the fact that pairwise 
additive potentials were employed and no adjustment of binary parameters to ternary data 
was carried out, all results of the present study on ternary systems are predictive. An 
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excellent agreement between the predictions and the experimental data was observed in 
most cases. Reliable predictions can also be expected for VLE of mixtures containing more 
than three components. 

Compared to the Peng-Robinson EOS, molecular modeling and simulation is found to 
yield superior predictions. 

Due to their numerical efficiency and accuracy, the molecular mixture models studied 
here are also well suited for simulations on a larger scale to investigate processes like 
evaporation, adsorption, flow etc. 
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SUPPORTING INFORMATION 

The supporting information contains the full numerical simulation data in comparison 
to the experiment for all studied ternary mixtures. Furthermore, the data is presented in 
ternary vapor-liquid phase diagrams. Note that for five systems figures are omitted, due to 
the fact that they either could not be simulated (N2 + R13 + R14) or only a single bubble 
point is available from experiment (R12 + R113 + R152a, R22 + R124 + R142b, R23 + 
R113 + R114 and R125 + R134a + R161). 



This material is available free of charge via the Internet at http://pubs.acs.org 
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APPENDIX, SIMULATION DETAILS 



The technical simulation details of the present calculations are similar to those published 
in [H [5]. A center-center cut-off radius of 17.5 A was used for the explicit evaluation of 
the intermolecular interactions. The Lennard- Jones tail corrections for internal energy, 
pressure, and chemical potential were calculated employing angle averaging as proposed 
by Lustig [Tlj. Long-range corrections for the dipolar part of the potential models were 
made with the reaction field method [73, [75]. The quadrupolar interaction needs no long 
range correction as it disappears by orientational averaging. The same holds for the mixed 
polar interaction between dipoles and quadrupoles, cf. Weingerl et al. [77] . 

VLE were obtained with the Grand Equilibrium method [78]. Depending on thermo- 
dynamic conditions, two levels of computational effort were employed: 

(A ) In simple cases (e.g. CH4 + CO2 + CaHg, CO2 + R142b + R152a and R13 + R14 + 
R23) VLE can be obtained with small statistical uncertainties sampling = 864 molecules 
for the liquid phase and about 500 molecules for the vapor phase. Liquid simulation runs 
were carried out using molecular dynamics with 200 000 time steps, vapor simulation runs 
were performed using the Monte Carlo technique with 200 000 cycles. Within one cycle, 
N attempts to translate or rotate, and two attempts to insert or delete molecules were 
sampled. The chemical potentials were calculated by Widom's insertion technique [71] 
using 3456 test molecules each time step. 

(B) In difficult cases (e.g. Ar + N2 + O2, RIO + R20 + R30 and R30 + R30B1 + 
R30B2), where experimental data is present only for highly dense strongly polar liquid 
phases and the vapor pressure is usually low, the more elaborate gradual insertion scheme 
had to be employed to obtain the chemical potentials in the liquid. 

The gradual insertion method is an expanded ensemble method [TH] based on the Monte 

11 



Carlo technique. The version as proposed by Nezbeda and Kolafa [5U], extended to the 
NpT ensemble [81], was used in case (B). In comparison to Widom's insertion technique, 
where full molecules are inserted into the fluid, gradual insertion introduces one fluctuating 
molecule that undergoes changes in a predeflned set of discrete states of coupling with all 
other molecules constituting the fluid. Preferential sampling is done in the vicinity of 
the fluctuating molecule. This concept leads to considerably improved accuracy of the 
residual chemical potential. Gradual insertion simulations were performed with N = 864 
molecules in the liquid phase. Starting from a face-centered cubic lattice arrangement, 
every simulation run was given 5000 Monte Carlo cycles to equilibrate. Data production 
was performed over 100 000 Monte Carlo cycles. One Monte Carlo cycle is deflned here 
as trial translations, (2/3) N trial rotations, and one trial volume change. Further 
simulation parameters for runs with the gradual insertion method were taken from Vrabec 
et al. [S]. 
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LIST OF SYMBOLS 
Latin Letters 



kij binary parameter of the Peng-Robinson equation of state 

L elongation 

p pressure 

Q quadrupolar momentum 

T temperature 

X mole fraction in liquid phase 

y mole fraction in vapor phase 



Greek Letters 



e Lennard- Jones energy parameter 

jj, dipolar momentum 

^ binary interaction parameter 

(7 Lennard- Jones size parameter 



Subscripts 



A related to component A 

B related to component B 

C related to component C 

i related to component i 

ij related to components i and j 

j related to component j 
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Superscripts 



exp experimental data 

sim simulation data 



Abbreviations 



ICLJ 


one-center Lennard- Jones 




ICLJD 


one-center Lennard- Jones plus 


point dipole 


2CLJ 


two-center Lennard- Jones 




2CLJD 


two-center Lennard- Jones plus 


point dipole 


2CLJQ 


two-center Lennard- Jones plus 


point quadrupole 


DDB 


Dortmunder Datenbank 




EOS 


equation of state 




VLE 


vapor-liquid equilibria 





14 



References 

[1] Vrabec, J.; StoU, J.; Hasse, H. J. Phys. Chem. B 2001, 105, 12126. 

[2] StoU, J.; Vrabec, J.; Hasse, H. J. Chem. Phys. 2003, 119, 11396. 

[3] Streett, W. B.; Tildesley, D. J. Proc. Roy. Soc. Lond. A 1977, 355, 239. 

[4] Vrabec, J.; Huang, Y.-L.; Hasse, H. Fluid Phase Equilib. 2009, in press. 

[5] StoU, J.; Vrabec, J.; Hasse, H. AIChE J. 2003, 49, 2187. 

[6] Vrabec, J.; StoU, J.; Hasse, H. Mol. Sim. 2005, 31, 215. 

[7] StoU, J. Molecular Models for the Prediction of Thermophysical Properties of Pure 
Fluids and Mixtures; Fortschritt-Berichte VDI; VDI-Verlag: Diisseldorf, 2005; Reihe 
3, Vol. 836. 

[8] Carrero-MantUla, J.; Llano- Restrepo, M. Fluid Phase Equilib. 2003, 208, 155. 

[9] Potoff, J. J.; Siepmann, J. I. AIChE J. 2001, 47, 1676. 
[10] Kamath, G.; Potoff, J. J. Fluid Phase Equilih. 2006, 246, 71. 
[11] Hansen, N.; Agbor, F. A. B.; Keil, F. J. Flmd Phase Equilib. 2006, 259, 180. 
[12] Liu, A.; Beck, T. L. J. Phys. Chem. B 1998, 102, 7627. 
[13] Nath, S. K.; Banaszak, B. J.; de Pablo, J. J. Macromolecules 2001, 34, 7841. 
[14] Lisal, M.; Smith, W. R.; Nezbeda, I. AIChE J. 2000, 46, 866. 
[15] Van't Hof, A. PhD thesis, TU Delft, 2005. 

15 



[16] Tsang, P. C; White, O. N.; Perigard, B. Y.; Vega, L. F.; Panagiotopoulos, A. Z. 
Flmd Phase Equilih. 1995, 107, 31. 

[17] de Miguel, E.; Telo da Gama, M. M. J. Chem. Phys. 1997, 107, 6366. 

[18] Escobedeo, F. A. J. Chem. Phys. 1999, 110, 11999. 

[19] Sadus, R. J. Fluid Phase Equilih. 1999, 157, 169. 

[20] Errington, J. R.; Shen, V. K. J. Chem. Phys. 2005, 123, 164103. 

[21] Attwood, B. C; Hall, C. K. AIChE J. 2008, 54, 1886. 

[22] Peng, D. Y.; Robinson, D. B. Ind. Eng. Chem. Fundam. 1976, 15, 59. 

[23] Dortmunder Datenbank, Mixture Properties, Version 1.3.0.211, 2004. 

[24] Miller, R. C; Kidnay, A. J.; Hiza, M. J. AIChE J. 1973, 19, 145. 

[25] Christiansen, L. J.; Fredenslund, Aa. Cryogenics 1974, I4, 10. 

[26] Elshayal, I. M.; Lu, B. C. Y. Cryogenics 1971, 11, 285. 

[27] Fastovsky, W. G.; Petrovsky Yu. W. Zh. Fiz. Khim. 1957, 31, 836. 

[28] Kremer, H.; Knapp, H. Fluid Phase Equilih. 1983, 11, 289. 

[29] Xu, N.; Dong, J.; Wang, Y.; Shi, J. Fluid Phase Equilih. 1992, 81, 175. 

[30] Trappehl, G.; Knapp, H. Cryogenics 1987, 27, 696. 

[31] Okambawa, R.; St-Arnaud, J. M.; Bose, T. K.; Le Noe, O. Fluid Phase Equilih. 
1997, 134, 225. 

[32] Knapp, H.; Yang, X.; Zhang, Z. Fluid Phase Equilih. 1990, 54, 1. 

16 



[33] Hsi, C; Lu, B. C. Y. Can. J. Chem. Eng. 1971, 49, 140. 

[34] Zenner, G. H.; Dana, L. I. Chem. Eng. Progr. Symp. Ser. 1963, 59, 36. 

[35] Zhang, Z.; Guo, L.; Yang, X.; Knapp, H. J. Chem. Ind. Eng. (China) 1999, 50, 392. 

[36] Al Wakeel, I. M. PhD thesis, Technische Universitat Berlin, 1976. 

[37] Wang, X.; Wang, Y.; Shi, J.; Lu, B. C. Y. J. Chem. Eng. Data 1991, 36, 436. 

[38] Zaharov, N. D.; Litvinov, A. S.; Grezin, A. K.; Gromov, E. A. Zh. Fiz. Khim. 1976, 
9, 162. 

[39] Bian, B. thesis, 1992. 

[40] McCurdy, J. L.; Katz, D. L. Ind. Eng. Chem. 1944, 36, 674. 

[41] Rulewicz, G.; Schuberth, H.; Leibnitz, E. J. Prakt. Chem. 1968, 37, 122. 

[42] Xie, R.; Liu, Y.; Yu, R; He, C; Qin, H.; Dong, B. Petrochem. Tech. 1981, 26, 144. 

[43] Chareton, A.; Valtz, A.; Laugier, S.; Richon, D.; Renon, H. J. Chem. Eng. Data 
1990, 35, 162. 

[44] Valtz, A.; Laugier, S.; Richon, D. J. Chem. Eng. Data 1987, 32, 397. 

[45] Proust, P. C.; Stein, F. P. Fluid Phase Equilib. 1979, 3, 313. 

[46] Bouchot, C.; Richon, D. Ind. Eng. Chem. Res. 1998, 37, 3295. 

[47] Cao, W.; Yu, H.; Wang, W. J. Chem. Ind. Eng. (China) 1997, 48, 136. 

[48] Laugier, S.; Richon, D.; Renon, H. Fluid Phase Equilib. 1994, 93, 297. 

[49] Apelblat, A.; Wisniak, J.; Tamir, A. J. Chem. Eng. Data 1981, 26, 144. 

17 



[50] Kobayashi, M.; Nishiumi, H. Fluid Phase Equilib. 1998, I44, 191. 

[51] Kruse, H.; Chen, J.; Koester, R.; Herres, G.; Buschmeiner, M.; Gorenflo, D. DKV- 
Tagungsberichte 1996, 23. Jahrgang, 67. 

[52] Kul, I.; DesMarteau, D. D.; Beyerlein, L. A. Fluid Phase Equilib. 2000, 173, 263. 

[53] Kang, Y. W.; Lee, Y. Y. J. Chem. Eng. Data 1996, 4I, 303. 

[54] Deiters, U. K. Fluid Phase Equilib. 1997, 132, 265. 

[55] StoU, J.; Vrabec, J.; Hasse, H.; Fischer, J. Fluid Phase Equilib. 2001, 179, 339. 

[56] StoU, J.; Vrabec, J.; Hasse, H. Fluid Phase Equilib. 2003, 209, 29. 

[57] Vrabec, J.; Kedia, G. K.; Hasse, H. Cryogenics 2005, 45, 253. 

[58] Vrabec, J.; Kumar, A.; Hasse, H. Fluid Phase Equilib. 2007, 258, 34. 

[59] Fernandez, G. A.; Vrabec, J.; Hasse, H. Int. J. Thermophys. 2004, 25, 175. 

[60] Fernandez, G. A.; Vrabec, J.; Hasse, H. Fluid Phase Equilib. 2004, 221, 157. 

[61] Fernandez, G. A.; Vrabec, J.; Hasse, H. Int. J. Thermophys. 2005, 26, 1389. 

[62] Fernandez, G. A.; Vrabec, J.; Hasse, H. Mol. Sim. 2005, 31, 787. 

[63] Fernandez, G. A.; Vrabec, J.; Hasse, H. Cryogenics 2006, 46-1 711. 

[64] Schnabel, T.; Vrabec, J.; Hasse, H. J. Mol. Liq. 2007, 135, 170. 

[65] Gray, C. G.; Gubbins, K. E. Theory of molecular fluids, Vol. 1, Fundamentals; 
Clarendon Press: Oxford, 1984. 

18 



[66] Allen, M. P.; Tildesley, D. J. Computer simulations of liquids; Clarendon Press: 
Oxford, 1987. 

[67] Smith, J. M.; Van Ness, H. C; Abbott, M. M. Introduction to chemical engineering 
thermodynamics, 5th Edition; McGraw-Hill: New York, 1996. 

[68] Merseburger Datenbank. 

[69] Case, F. H.; Brennan, J.; Chaka, A.; Dobbs, K. D.; Friend, D. C; Frurip, D.; Gordon, 
P. A.; Moore, J.; Mountain, R. D.; Olson, J.; Ross, R. B.; Schiller, M.; Shen, V. K. 
Fluid Phase Equilib. 2007, 260, 153. 

[70] Andersen, H. C. J. Chem. Phys. 1980, 72, 2384. 

[71] Widom, B. J. Chem. Phys. 1963, 39, 2808. 

[72] Heyes, D. M. Mol. Sim. 1992, 8, 227. 

[73] Flyvbjerg, H.; Petersen, H. G. J. Chem. Phys. 1989, 91, 461. 

[74] Lustig, R. Mol. Phys. 1988, 65, 175. 

[75] Barker, J. A.; Watts, R. O. Mol. Phys. 1973, 26, 789. 

[76] Saager, B.; Fischer, J.; Neumann, M. Mol. Sim. 1991, 6, 27. 

[77] Weingerl, U.; Fischer, J. Fluid Phase Equilib. 2002, 202, 49. 

[78] Vrabec, J.; Basse, H. Mol. Phys. 2002, 100, 3375. 

[79] Shevkunov, S. V.; Martinovski, A. A.; Vorontsov-Velyaminov, P. N. High Temp. 
Phys. (USSR) 1988, 26, 246. 

[80] Nezbeda, I.; Kolafa, J. Mol. Sim. 1991, 5, 391. 

19 



[81] Vrabec, J.; Kettler, M.; Hasse, H. Chem. Phys. Lett 2002, 356, 431. 



20 



Table 1: List of the 35 components studied in the present work. The model parameters 
were taken from [H [2] . 



Component 


PAS 


Component 


PAS RN 


Non-polar, ICLJ 




R141b 


fCH3-CFCl2) 


1717-00-6 


Ar 




13965-95-2 


R142b 


CH3-CF2CI) 


75-68-3 


CH4 




74-82-8 


R143a 


CH3-CF3) 


420-46-2 


Dipolar, ICLJD 




R152a 


^CH3-CHF2) 


( 0-0 i -u 


R30 (CH2CI2) 


75-09-2 


R161 (CH2F-CH3) 


353-36-3 


R30B2 (CHsBra) 


74-95-3 


Quadrupolar, 2CLJQ 




R32 (CH2F2) 


75-10-5 


N2 




7727-37-9 


Dipolar, 2CLJD 




O2 




7782-44-7 


CO 




630-08-0 


CO2 




124-38-9 


Rll 


CFCI3) 


75-69-4 


C2H2 




74-86-2 


R12 


CF2CI2) 


75-71-8 


C2H4 




74-85-1 


R13 


CF3CI) 


75-72-9 


C2H6 




74-84-0 


R20 


'CHCI3) 


67-66-3 


RIO (CCI4) 


56-23-5 


R22 


'CHF2CI) 


75-45-6 


R14 (CF4) 


75-73-0 


R23 


'CHF3) 


75-46-7 


R113 (CFCI2-CF2CI) 


76-13-1 


R30B1 (CH2BrCl) 


74-97-5 


R114 (CF2CI-CF2CI) 


76-14-2 


R124 


(CHFCI-CF3) 


2837-89-0 


R150B2 (CH2Br-CH2Br) 


106-93-4 


R125 


CHF2-CF3) 


354-33-6 


RlllO 


;c2ci4) 


127-18-4 


R134a (CH2F-CF3) 


811-97-2 


R1120 


|CHCl=CCl2) 


79-01-6 


R140a (CCI3-CH3) 


71-55-6 
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Table 2: List of the 33 studied ternary mixtures and reference to literature on experimental 
VLE. For systems indicated with f, only bubble line data is available from experiment. 



Ar -f 


- CH4 + N2 


[21 


N2 + CO2 + C2H6 




R13 + R14 + R23 |15] 


Ar 4- 


- CH4 + CO 




N2 + CO2 + R12 


[36] 


R22 + R23 + R114 f |13] 


Ar 4- 


- CH4 + C2H6 


[26| 


N2 + CO2 + R22 


[37J 


R22 + R124 + R142b f |46] 


Ar + N2 + O2 


m 


N2 + R13 + R14 t 


[38J 


R22 + R142b + R152a |17| 


CH4 


+ N2 + CO 


[2H| 


CO2 + R22 + R142b 


[39J 


R23 + R113 + R114 t [IE] 


CH4 


+ N2 + C02 


m 


CO2 + R142b + R152a 


[39J 


R30 + R30B1 + R30B2 ^ 


CH4 


+ N2 + C2H6 


m 


C2H2 + C2H4 +C2H6 


[40j 


R32 + R125 + R134a [SO] 


CH4 


+ CO + CO2 


m\ 


RIO + R20 + R30 


[41J 


R32 + R125 + R143a [SI] 


CH4 


+ C02 + C2H6 




RIO + RlllO + R1120 


[42J 


R125 + R134a + R143a [S2] 


CH4 


+ C2H4 + C2He 


[33J 


Rll + R22 + R23 f 


[43J 


R125 + R134a + R161 f 


N2 + O2 + CO2 


[34J 


R12 + R113 + R152a f 


[44J 


R140a + R141b + R142b [M] 
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